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ABSTRACT 

We study the evolution of the cosmic star formation in the universe by computing 
the luminosity density (in the UV, B, J, and K bands) and the stellar mass density 
of galaxies in two reference models of galaxy evolution: the pure-luminosity evolution 
(PLE) model developed by Calura & Matteucci (2003) and the semi-analytical model 
(SAM) of hierarchical galaxy formation by Menci et al. (2002). The former includes 
a detailed description of the chemical evolution of galaxies of different morphological 
types; it does not include any number evolution of galaxies whose number density 
is normalized to the observed local value. On the other hand, the SAM includes a 
strong density evolution following the formation and the merging histories of the DM 
haloes hosting the galaxies, as predicted by the hierarchical clustering scenario, but 
it does not contain morphological classification nor chemical evolution. Our results 
suggest that at low-intermediate redshifts {z < 1.5) both models are consistent with 
the available data on the luminosity density of galaxies in all the considered bands. At 
high redshift the luminosity densities predicted in the PLE model show a peak due to 
the formation of ellipticals, whereas in the hierarchical picture a gradual decrease of the 
star formation and of the luminosity densities is predicted for z > 2.5. At such redshifts 
the PLE predictions tend to overestimate the present data in the B band whereas the 
SAM tends to underestimate the observed UV luminosity density. As for the stellar 
mass density, the PLE picture predicts that nearly 50% and 85% of the present stellar 
mass are in place at z ^ A and z ^ 1, respectively. According to the hierarchical 
SAM, 50% and 60% of the present stellar mass are completed at z ~ 1.2 and z = 1, 
respectively. Both predictions fit the observed stellar mass density evolution up to 
z = 1. At z > 1, the PLE and SAM models tend to overestimate and underestimate 
the observed values, respectively. We discuss the origin of the similarities and of the 
discrepancies between the two models, and the role of observational uncertainties (such 
as dust extinction) in comparing models with observations. 

Key words: Galaxies: formation and evolution; Galaxies: fundamental parameters. 



1 INTRODUCTION 

In the past few years a great deal of work appeared on the 
subject of galaxy formation and evolution. With the word 
"formation" usually one means the assembly of the bulk of 
the material (say > 50%) of the luminous part of a galaxy, 
namely the stars and the gas, within a sphere of radius of 
^30 kpc (Peebles 2003) . A reliable picture of galaxy forma- 
tion must be able to reproduce, at the same time, all (or as 
much as possible) of the available constraints, including col- 
ors and chemical abundances. Currently, the most intriguing 
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debate on galaxy evolution concerns how the formation of 
ellipticals and bulges occurred in the universe. In fact, the 
two main competing scenarios of galaxy evolution propose 
rather different conditions for the formation of spheroids. In 
the first scenario, ellipticals and bulges formed at high red- 
shift (e.g. 2 > 2 — 3) as the result of a violent burst of star 
formation following a "monolithic collapse" (MC) of a gas 
cloud. After the main burst of star formation, the galaxy lost 
the residual gas by means of a galactic wind and it evolved 
passively since then (Larson 1974, van Albada 1982, Sandage 
1986, Matteucci & Tornambe 1987, Arimoto & Yoshii 1987, 
Matteucci 1994). The monolithic collapse view, or better, 
the idea that spheroids formed quickly and at high red- 
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shift, is supported by a large set of observational evidences. 
Among them, of particular importance are the thinness of 
the Fundamental Plane (Djorgovski & Davis 1987, Renzini 
& Ciotti 1993, Bernardi et al. 1998, Kochanek et al. 2000, 
van Dokkum et al. 2001, Rusin et al. 2003, van Dokkum 
& Ellis 2003), the overabundance of Mg relative to Fe ob- 
served in the stars as well as the increase of the [Mg/Fc] 
ratio with galaxy luminosity (Pipino & Matteucci 2004 and 
references therein) , the tightness of the color-central veloc- 
ity dispersion and color-magnitude relation (Bower, Lucey 
& Ellis 1992, Kodama et al. 1999) observed for both cluster 
and field spheroids at high and low redshift, as well as the 
constancy of the number density of both spheroids and large 
discs observed up to z ~ 1 (Im et al. 1996, Lilly et al. 1998, 
Schade et al. 1999, Im et al. 2002). 

On the other hand, the hierarchical clustering (HC) picture 
is based on the Press & Schechter (1974) structure forma- 
tion theory, which has been developed mainly to study the 
behaviour of the dark matter. According to this theory, in 
a A-Cold dark Matter (ACDM)-dominated universe, small 
DM halos are the first to collapse, then interact and merge 
to form larger halos. The most uncertain assumption in the 
HC scenario concerns the behaviour of the baryonic matter, 
which is assumed to follow the DM in all the interaction and 
merging processes. In this framework, massive spheroids are 
formed from several merging episodes among gas-rich galax- 
ies, such as discs, occurring throughout the whole Hubble 
time. These mergers produce moderate star formation rates 
(SFRs), with massive galaxies reaching their final masses at 
more recent epochs than less massive ones (z < 1.5, White & 
Rees 1978, Kauffmann, White & Guiderdoni 1993, Baugh et 
al. 1998, Cole et al. 2000, Somerville et al. 2001, Menci et al. 
2002). The observational evidence in favor of the hierarchical 
galaxy formation is the apparent paucity of giant galaxies 
at high redshift (z ~ 1) as claimed by some authors (Barger 
et al. 1999, Kauffmann, Chariot & White 1996, Zepf 1997), 
the blue colors of some spheroids at low redshift, possibly 
ascribed to residual star formation activity induced by merg- 
ers (Franceschini et al. 1998, Menanteau et al. 1999), as well 
as the observations showing evidence for mergers in distant 
field and cluster galaxies (Bundy et al. 2004, van Dokkum 
et al. 2000) and the increase of the measured merging rate 
with redshift (Patton et al. 1997, Le Fevre et al. 2000, Con- 
selice et al. 2003). 

Recently, Calura & Matteucci (2003, hereinafter CM03) 
have developed a series of detailed chemical and spectro- 
photometric models for elliptical, spiral and irregular galax- 
ies, used to study the evolution of the luminous matter in 
the universe and the contributions that galaxies of different 
morphological types bring to the overall cosmic star forma- 
tion. It is worth noting that all these models reproduce the 
chemical abundances and abundance patterns in the afore- 
mentioned galaxies. 

In their scenario of pure-luminosity evolution (PLE), 
only the galaxy luminosities evolve, whereas the number 
densities are assumed to be constant and equal to the values 
indicated by the local B-band luminosity function (LF), as 
observed by Marzke et al. (1998). In this paper, we compare 
the cosmic star formation history as predicted by the PLE 
model of CM03 with the predictions of the hierarchical semi- 
analytic model (SAM) developed by Menci et al. (2002). We 
want to stress that the PLE model and the SAM do not 



represent the only alternatives to study galaxy evolution. 
For instance, several groups study the evolution of the cos- 
mic star formation by means of large-scale hydrodynamical 
simulations (e.g. Sringel & Hernquist 2003, Nagamine et al. 
2004), which are generally based on the ACDM cosmologi- 
cal model. However, representing the PLE and Menci SAM 
considered in this work two rather opposite scenarios and 
providing rather extreme predictions, they may be helpful 
to constrain the parameter space also for other galaxy for- 
mation models. Furthermore, we want to stress that not nec- 
essarily the two scenarios are in contradiction, since the HC 
was devised for the DM whereas the PLE for the baryonic 
matter. As some observational evidence seems to indicate, it 
is in fact possible that, although DM halo formation is hier- 
archical, the baryonic matter evolved in an anti-hierarchical 
fashion, in the sense that larger galaxies are older than small 
ones (Matteucci 1994, Pipino & Matteucci 2004). By com- 
paring the model predictions with a large set of observational 
data, we aim at inferring whether the two main competing 
scenarios can be disentangled on the basis of the current ob- 
servations. The novelty with respect to the paper by CM03 is 
the incorporation of dust extinction in the PLE model, with 
important consequences on the predicted behaviour of the 
luminosity of galaxies at short wavelengths, i.e. in the UV 
and B photometric bands. This paper is organized as follows: 
in sections 2 and 3, we describe the pure-luminosity evolu- 
tion model as developed by CM03 and the SAM by Menci et 
al. (2002), respectively. In section 4 we present our results, 
and in section 5 we draw the conclusions. Unless otherwise 
stated, throughout the paper we use a ACDM cosmological 
model characterized by f2o = 0.3, f^A = 0.7 and h — 0.65. 



2 THE CMOS PURE-LUMINOSITY 
EVOLUTION MODEL 

The PLE models developed by CM03 consist of chemi- 
cal evolution models for galaxies of different morphologi- 
cal types (ellipticals, spirals, irregulars), used to calculate 
metal abundances and star formation rates (SFRs), and by 
a spectro-photometric code used to calculate galax;y spectra, 
colors and magnitudes by taking into account the chemi- 
cal evolution. Detailed descriptions of the chemical evolu- 
tion models for galaxies of different morphological types can 
be found in Matteucci & Tornambe (1987) and Matteucci 
(1994) for elliptical galaxies, Chiappini et al. (1997, 2001) 
for the spirals and Bradamante et al. (1998) for irregular 
galaxies. We assume that the category of galactic bulges is 
naturally included in the one of elliptical galaxies. Our as- 
sumption is supported by the similar features characterizing 
bulges and ellipticals: for instance, both are dominated by 
old stellar populations and respect the same fundamental 
plane (Binney & Merrifield 1998, Renzini 1999). This indi- 
cates that they are likely to have a common origin, i.e. both 
are likely to have formed on very short timescales and a long 
time ago, and we will refer to both ellipticals and bulges as 
to the "spheroids" . 

In our picture, spheroids form as a result of the rapid col- 
lapse of a homogeneous sphere of primordial gas where star 
formation is taking place at the same time as the collapse 
proceeds. Star formation is assumed to halt as the energy 
of the ISM, heated by stellar winds and SN explosions, bal- 
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ances the binding energy of the gas. At this time a galactic 
wind occurs, sweeping away almost all of the residual gas. 
By means of the galactic wind, ellipticals enrich the inter- 
galactic medium (IGM) with metals. 

For spiral galaxies, the adopted model is calibrated in order 
to reproduce a large set of observational constraints for the 
Milky Way galaxy (Chiappini et al. 2001). The Galactic disc 
is approximated by several independent rings, 2 kpc wide, 
without exchange of matter between them. In our picture, 
spiral galaxies are assumed to form as a result of two main 
infall episodes. During the first episode, the halo and the 
thick disc are formed. During the second episode, a slower 
infall of external gas forms the thin disc with the gas accu- 
mulating faster in the inner than in the outer region ("inside- 
out" scenario, Matteucci & Frangois 1989). The process of 
disc formation is much longer than the halo and bulge for- 
mation, with time scales varying from ~ 2 Gyr in the inner 
disc to ~ 8 Gyr in the solar region and up to 10 — 15 Gyr in 
the outer disc. 

In this case, at variance with Chiappini et al. (2001) CMOS 
assume a Salpeter (1955) IMF, instead of the Scalo (1986) 
IMF. This choice is motivated by the fact that a Scalo or a 
Salpeter IMF in spirals produce very similar results in the 
study of the luminosity density evolution, and also by the 
fact that we aim to test the hypothesis of a universal IMF 
(see also Calura & Matteucci 2004). Another difference be- 
tween the Chiappini et al. (2001) model and ours concerns 
the elimination of the star formation threshold, motivated 
by the fact that its effects are appreciable only on small 
scales, i.e. in the chemical evolution of the solar vicinity and 
of small galactic regions, whereas our aim is to study star 
formation in galactic discs on global scales. Finally, irregular 
dwarf galaxies are assumed to assemble from continuous in- 
fall of gas of primordial chemical composition, until masses 
in the range ~ 10* — 6 x 1O^M0 are accumulated, and to 
produce stars at a lower rate than spirals. 
Let Gi be the fractional mass of the element i in the gas 
within a galaxy, its temporal evolution is described by the 
basic equation: 



xpit) = vGit) 



(2) 



G, 



-i>{t)X,{t) + R,(t) + (GO.n/ - (Gi)out 
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where Giii) = Ug{t)Xi{t)/atot is the gas mass in the form 
of an element i normalized to a total initial mass Mtot- 
The quantity Xi{t) = Gi{t)/G{t) represents the abundance 
in mass of an element i, with the summation over all ele- 
ments in the gas mixture being equal to unity. The quantity 
G{t) — ag{t)/atot is the fractional mass of gas present in 
the galaxy at time t. ipit) is the instantaneous star forma- 
tion rate (SFR), namely the fractional amount of gas turning 
into stars per unit time; Ri{t) represents the returned frac- 
tion of matter in the form of an element i that the stars 
eject into the ISM through stellar winds and SN explosions; 
this term contains all the prescriptions regarding the stellar 
yields and the SN progenitor models. The two terms (Gi)i„/ 
and {Gi)out account for the infalling external gas from the 
IGM and for the outflow, occurring by means of SN driven 
galactic winds, respectively. The main feature characteriz- 
ing a particular morphological galactic type is represented 
by the prescription adopted for the star formation history. 
In the case of elliptical and irregular galaxies the SFR ip^t) 
(in Gyr~^) has a simple form and is given by: 



The quantity v is the efficiency of star formation, 
namely the inverse of the typical time scale for star forma- 
tion and for ellipticals and bulges is assumed to be ~ 10— 15 
Gyr~^ (Matteucci 1994). In the case of spheroids, v is as- 
sumed to drop to zero at the onset of a galactic wind, which 
develops as the thermal energy of the gas heated by super- 
nova explosions exceeds the binding energy of the gas (Ari- 
moto & Yoshn 1987, Matteucci & Tornambe 1987). This 
quantity is strongly influenced by assumptions concerning 
the presence and distribution of dark matter (Matteucci 
1992); for the model adopted here a diffuse {Re/Rd=0.1, 
where R^ is the effective radius of the galaxy and Rd is the 
radius of the dark matter core) but massive {Mdark/MLum = 
10) dark halo has been assumed. 

In the case of irregular galaxies we have assumed a con- 
tinuous star formation rate always expressed as in (2), but 
characterized by an efficiency lower than the one adopted 
for ellipticals, i.e. v — 0.01 Gyr~^ 
In the case of spiral galaxies, the SFR expression is: 



tl^{r,t) = uat^t{r,t)ag'^{r,t) 



(3) 



where fei = 0.5 and k2 = 1.5 (see Matteucci & Frangois 
1989, Chiappini et al. 1997). For massive stars (M > 8Mq) 
we adopt nucleosynthesis prescriptions by Nomoto et al. 
(1997a), the yields by van den Hoeck & Groenewegen (1997) 
for low and intermediate mass stars (0.8 < M/Mq < 8) and 
those of Nomoto et al. (1997b) for type I a SNe. 
For all galaxies, we assume a Salpeter IMF, expressed by 
the formula: 



(j)om' 



-(l + x) 



(4) 



with X — 1.35, being the mass range 0.1 < m/rriQ < 100. To 
calculate galaxy colors and magnitudes, we use the photo- 
metric code by Bruzual & Chariot (2003, hereinafter EC). 
However, we have implemented the EC code by taking 
into account the evolution of metallicity in galaxies (Calura 
2004). Dust extinction is also properly taken into account. 
The chosen geometrical dust distribution plays an impor- 
tant role in the modelling of dust attenuation in galaxies: 
usually, the "screen" and "slab" dust distributions represent 
the two most extreme cases. In the screen model, the dust 
is distributed along the line of sight of the stars, whereas in 
the slab model the dust has the same distibution as stars. 
The main difference between the screen and slab dust dis- 
tributions is the expression of the attenuation factor, which 
in the former case is given by: 

ascreen = OXp (-r(A)) (5) 

whereas in the latter case it is given by: 

asia6 = [l-oxp(-r(A))]/r(A) (6) 

(Totani & Yoshii 2000), where t(A) is the optical depth of 
the dust. In this case, we adopt the "screen" geometric dis- 
tribution which, according to UV and optical observations 
of local starburst galaxies, is to be considered favored over 
the "slab" model (Calzetti et al. 1994). The absorbed flux 
7a (A) of a stellar population behind a screen of dust is given 
by: 



J4A)=77(A)exp(-r(A)) 



(7) 
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(Calzetti 2001), where //(A) represents the intrinsic, unob- 
scured flux at the wavelength A. 

We assume that the optical depth is proportional to the 
column density N{g) and to the metallicity Z of the gas, 
according to: 



-{\)^Ck{\)N{g)Z 



(8) 



where fc(A) is the extinction curve. For spiral galaxies, we 
adopt the extinction curve derived by Seaton (1979) for the 
Milky Way (MW) galaxy. Such a choice is motivated by 
the fact that we assume that, as far as the chemical and 
photometric features are concerned, the Milky Way Galaxy 
represents an average spiral. Local starburst galaxies are 
generally characterized by extinction curves slightly differ- 
ent from the ones of the MW (Calzetti 1997, 2001) and are 
better modelled by the expression found by Calzetti (1997). 
We assume that in the starbursts occurring in elliptical and 
irregular galaxies the dust follows an attenuation law similar 
to the one estimated by Calzetti (1997) for local starbursts. 
The constant C in equation (6) is chosen in order to repro- 
duce the Milky Way average V-band extinction of Ay = 0.17 
(Schlegel, Finkbeiner & Davis 1998). 

The galaxy densities of the various morphological types are 
normalized according to the local B-band luminosity func- 
tion observed by Marzke et al. (1998). A scenario of pure 
luminosity evolution has been assumed, namely that galax- 
ies evolve only in luminosity and not in number. This is 
equivalent to assume that the effects of galaxy interactions 
and mergers are negligible at any redshift. Such a picture 
can account for many observables, such as the evolution of 
the galaxy luminosity density in various bands and the cos- 
mic supernova rates (CMOS). At redshift larger than zero 
the absolute magnitudes are calculated according to: 



Mb{z) = Mb{z = Q)+ 2.5log{ 






(9) 



where Mb{z = 0) and Mb{z) wee the absolute blue magni- 
tudes at redshift and z, respectively, E\(z) dX is the energy 
per unit time radiated at the rest-frame wavelength A by the 
galaxy at redshift z, and Rb{X) is the response function of 
the rest-frame B band. The second term on the right side 
of equation 7 represents the evolutionary correction (EC), 
i.e. the difference in absolute magnitude measured in the 
rest frame of the galaxy at the wavelength of emission (Pog- 
gianti 1997). 
For the LF, we assume a Schechter (1976) form, given by: 

cE'(M) dM = 0.4 /n(10) $* e'^ X"+^ dM (10) 

where X = L/L* = iqO-Mm-m')^ ^. ^^,-^ j^ ^j^^ character- 
istic magnitude (luminosity) and is a function of redshift, 
whereas $* and a are the normalization and the faint-end 
slope, respectively, and are assumed to be constant. 
In bands other than B we assume that the LF shape is the 
same as in the B band and we calculate the LF in the given 
band (X) transforming the absolute magnitudes according 
to the rest-frame galaxy colors as predicted by the spec- 
trophotometric model: 



Mx = Mb + {X - B] 



rf 



(11) 



The LD per unit frequency in a given band (centered at the 
wavelength A) and for the fc— morphological type is: 



p\,k= l ^k{,Lx){Lx/Ll)dLx 



(12) 



The total LD is given by the sum of the single contributions 
of spheroids, spirals and irregulars. 

The stellar mass densities for galaxies of the k~th morpho- 
logical type are p*,fe and are calculated as: 



P*,k = pB,k ■ {M,/L)B,k 



(13) 



where pB,k is the predicted B luminosity density, whereas 
{Mt/L)B,k is the predicted stellar mass to light ratio for 
the A:— th galactic morphological type. All the galaxies are 
assumed to start forming stars at the same redshift Zf — 5. 



3 THE SAM MODEL 

In semianalytical models the galaxy mass distribution is de- 
rived from the merging histories of the host DM haloes, un- 
der the assumption that the galaxies contained in each halo 
coalesce into a central dominant galaxy if their dynamical 
friction timescale is shorter than the halo survival time. The 
surviving galaxies (commonly referred to as satellite galax- 
ies) retain their identity and continue to orbit within the 
halo. The histories of the DM condensations rely on a well 
established framework (the extended Press & Schechter the- 
ory, EPST, see Bower 1991; Bond et al. 1991; Lacey & Cole 
1993). However, the recipe concerning the galaxy fate inside 
the DM haloes is guided by a posteriori consistency with the 
outputs of high-resolution N-body simulations. The SAM in- 
cludes the main dynamical processes taking place inside the 
host DM halos, namely dynamical friction and binary ag- 
gregations of satellite galaxies. The evolution of the galaxy 
mass distribution is calculated by solving numerically a set 
of evolutionary equations (Poll et al. 1999). 
The link between stellar evolution and the dynamics fol- 
lows a procedure widely used in semianalythic models. The 
baryonic content (ilb/fim) rn of the galaxy is divided into 
(1) a hot phase with mass nih at the virial temperature 
T — (1/2) /im// w^/fc (niH is the proton mass and p is the 
mean molecular weight), (2) into a cold phase with mass rric 
able to radiatively cool within the galaxy survival time, and 
the stars (3) (with total mass m*) forming from the cold 
phase on a time scale r, . Initially, all baryons are assigned 
to the hot phase. 

Also in this case, we compute galaxy spectra and luminosi- 
ties by means of the spectrophotometric code developed by 
Bruzual & Chariot (2003). The integrated stellar emission 
Sx{v, i) at the wavelength A for a galaxy of circular velocity 
V at the time t is computed by convolving with the spectral 
energy distribution <j)x obtained from population synthesis 
models: 



Sx{v,t)= dt' <j)x{t-t')rn,{v,t') 
Jo 



(14) 



(px is taken from Bruzual & Chariot (2003), with a Salpeter 
IMF. The metallicity is calculated by assuming a constant 
effective yield. The average galaxy metallicity varies between 
Z ~ 0.003 and Z ~ 0.01, in agreement with results of other 
SAMs (e.g.. Cole et al. 2000). To calculate galactic spectra, 
we use simple stellar populations (SSPs) at fixed metallicity 
Z = 0.004. The use of the SSPs at Z = 0.008 would produce 
very small variations in our results, certainly of negligible 
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entity with respect to the observational errors. 
The dust extinction affecting the above luminosities is com- 
puted assuming the dust optical depth to be proportional 
to the metallicity Z^oid of the cold phase and to the disk 
surface density, so that for the V- band rv oc rric Zcoid/i^r^- 
The proportionality constant is taken as a free parameter 
chosen to fit the bright end of the local LF. This fact yields, 
for the proportionality constant, the value 3.5 Mq^ pc^ with 
the stellar yield producing a solar metallicity for & v — 220 
km/s galaxy. Physically, this recipe for computing dust ex- 
tinction is identical to the one used for the PLE model (eq. 
6). To compute the extinction in the other bands, we use 
the extinction law of Calzotti (1997). 



4 PLE VS SAMS : RESULTS 

4.1 The SFR density 

In Figure 1, we show the evolution of the cosmic SFR den- 
sity as a function of redshift as predicted in the framework of 
the two scenarios. The two curves have very different shapes: 
the PLE scenario predicts a peak at redshift « ~ 5 due to 
starbursts in spheroids (CMOS), followed by a flat behavior 
between z ~ 4.2 and z ~ 3 due to star formation in spiral 
galaxies. The maximum SF in spirals cause a smaller peak 
of star formation at z = 2, and these galaxies are the re- 
sponsible for the decline of the SFR density between z — 2 
and z = 0. 

The hierarchical SAM model by Menci et al. (2002) produces 
a curve characterized by a weak increase between z — 5 and 
2^3, then becomes constant between z ~ 3 and z ~ 2 and 
finally starts to decrease at z < 2 down to z = 0. Between 
z = 2 and z = 0, the SAM model predicts a higher amount 
of SF than the PLE one. 



4.2 The galaxy luminosity density 

In Figure 2, we show the redshift evolution of the luminos- 
ity density in the rest-frame K (lower panel) and J (upper 
panel) bands, as predicted by the PLE (solid lines) and by 
the SAM (dashed lines) , compared to a set of observational 
data by various authors. 

The K band, centered at A = 2.2/i, is dominated by long 
-lived, low mass stars. The light emitted in this band is 
unaffected by dust extinction. At z > 2, the two curves 
have dramatically different behaviours: the PLE shows the 
peak due to ellipticals, whereas the SAM curve has a broad 
peak centered at z ~ 2. On the other hand, it is compelling 
how similar the curves are at z < 2. At z < 1 we show 
the observational data by Pozzetti et al. (2003) and Cohen 
(2002), in substantial agreement with one another. In this 
redshift range, both curves show broadly a good agreement 
with the observational data. The PLE scenario predicts a 
slightly higher LD at z = 0, mainly due to the higher num- 
ber of old stars (hence to redder galaxy colors) than the 
hierarchical picture. From the current set of observational 
data in the K band, it is practically impossible to distin- 
guish between the two opposite galaxy formation scenarios. 
Rest-frame Near Infrared deep galaxy surveys aimed at de- 
tecting faint sources, possibly located at high redshift, could 
provide us with fundamental hints to disentangle between 



the PLE and the hierarchical scenario. In fact, if there were 
an epoch when the bulk of spheroidal galaxies is forming, 
the K-band LD would show a peak centered at the redshift 
corresponding to that epoch. On the other hand, if massive 
galaxy formation is distributed throughout an extended pe- 
riod, no peak in the K band LD should be visible at high 
redshift. These results indicate that the study of the evolu- 
tion of the K band luminosity density at redshift larger than 
2 could represent the most direct observational strategy to 
establish the best scenario of galaxy formation. 

Similar conclusions can be drawn in the J band, domi- 
nated both by relatively old stars experiencing the red giant 
branch phase and by young main -sequence stars and in very 
similar fractions (Bruzual 2003). 

The above results, concerning the luminosity density in 
bands where the contribution of long-lived stars is relevant, 
show that the PLE and SAM models correctly predict the 
total amount of stars formed by z ~ 0, a conclusion con- 
firmed by our analysis of the stellar mass density (see be- 
low, sect. 4.3). The difference between the two scenarios is 
related to the rate of star formation during the cosmic time, 
which is better probed in the UV and B bands, where the 
contribution from massive, young stars is dominant. 

Figure 3 shows the evolution of the rest-frame UV and B 
luminosity density, as predicted by the PLE (solid curves) 
and SAM models (dashed curves). In this case, the theo- 
retical LDs have been calculated at 1400A and have been 
compared with data measured at various wavelengths, rang- 
ing from 1500 to 2800 A (see caption to Fig.3 for further 
details). In the two upper panels, the theoretical predictions 
are not corrected for dust extinction, whereas in the two 
lower panels the curves take into account also corrections for 
dust-extinction. Looking at the upper left panel it is possible 
to see how, once dust correction is not taken into account, 
in the UV band the PLE scenario predicts a strong peak at 
redshift 5. This peak is due to star formation in spheroids, 
which is absent in the hierachical scenario of Menci et al. 
(2002). On the other hand, the SAM curve shows a broad 
peak, centered at redshift ~ 2.5. Another difference concerns 
the predicted evolution at redshift < 1, where the curve from 
the SAM is constantly higher than the PLE one. This re- 
flects the fact that the SAM model predicts a higher amount 
of star formation occurring at z < 1 than the PLE curve; 
this is mainly due to the contribution of small- mass galaxies, 
which retain a relevant fraction of their gas down to small z, 
while the massive galaxy population, originated from clumps 
formed at high z in high-density regions, has already con- 
sumed most of the available cold gas reservoir. 
The curves calculated in the B band (upper right panel) 
show a behaviour very similar to the UV band, since both 
are dominated by the same types of stars, i.e. the youngest 
and the most massive ones. Both bands are sensitive to dust 
extinction, but in a different way: a comparison with the ob- 
servations can be discussed only after having corrected the 
curves for dust obscuration. 

In the lower left panel of Figure 3, the predicted UV lumi- 
nosity densities have been corrected for dust extinction. A 
very important result regarding the UV luminosity density 
predicted by the PLE scenario is that, once dust effects are 
properly taken into account, the peak at z ~ 5 due to el- 
lipticals appears considerably reduced, with the PLE curve 
showing a flat behaviour as the observational data by Pas- 
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Figure 1. The global SFR density versus redshift as predicted by the PLE model (solid line) and by the hierarchical SAM (dotted line) 
of Menci et al. 2002. 
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Figure 2. Luminosity density evolution in the rest-frame J (upper panel) and K (lower panel) bands as predicted by the PLE model by 
CMOS (solid curves) and by the hierarchical SAM of galaxy formation by Menci et al. (2002,das/ie(i curves), and as observed by Lilly et 
al. (1996, open circles), Pozzetti et al. (2003, solid triangles), Gardner et al. (1997, three-tips stars), Cohen (2002, solid squares). 
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Figure 3. Luminosity density evolution in the rest-frame UV and B bands as predicted by our PLE models {solid curves) and by the 
hierarchical SAM of galaxy formation by Menci et al. (2002) {dashed curves) and as observed by varour authors. For the UV band, the 
theoretical curves have been calculated at rest-frame 1400 A. UV band observations: Cowie et al. (1999, 2500A, four-tips stars), Pascarelle 
et al. (1998, 1500 A, solid hexagons), Steidel et al. (1999, 1500 A, open hexagons), Treyer et al. (1998, 2000 A, cross), Massarotti et al. 
(2001, 1500 A, five-tips stars), Giavalisco et al. (2004, 1500 A solid pentagons), Lilly et al. (1996, 2800 A, open circles), Connolly et al. 
(1997, 2800 A, open squares), Lanzetta et al. (2002, 1500 A, solid diamonds, plotted for different values of the parameters involved in 
their measure). B band observations: Ellis et al. (1996, 4400 A, open triangles), Dickinson et al. (2003, 4500 A, solid circles), Rudnick 
et al. (2003, 4400 A, open diamonds), Connolly et al. (1997, 4400 A, open squares), Lilly et al. (1996, 4400 A, open circles). Wolf et al. 
(2003, 4560 A, open pentagons). In the two upper panels, the theoretical curves are not corrected for dust extinction. In the two lower 
panels, the curves take into account dust extinction corrections. 



carelle et al. (1998) and Steidel et al. (1999). This means 
that, as suggested by CMOS, if the bulk of the star forma- 
tion in the high-redshift universe occurred in sites highly ob- 
scured by dust, most of it would be invisible for rest-frame 
UV surveys (see also Franx et al. 2003). Of great interest 
would be the study of the IR/submm luminosity density, 
which would be considerably enhanced by the re-emission 
by dust of all the UV absorbed flux, and which is deferred 
to a forthcoming paper. It is also important to note that at 
redshift > 4, the dust-corrected prediction from the hierar- 
chical model is critical: at very high redshift, the unobscured 
UV luminosity density (and hence the amount of star for- 
mation) is probably underestimated by the SAM by a factor 
of 3 or more, although the scatter in the data is too large 
to draw firm conclusions. However, recent independent anal- 
ysis (Fontana et al. 2003b, Menci et al. 2004) have shown 
that when only the bright galaxy population is selected, the 
paucity of the predicted UV luminosity density compared 
with observations is more clearly revealed, confirming that 
at those z some fundamental process must be at work, such 
as bursts of star formation with a rate higher than that 
predicted by standard SAMs. Such a process could be con- 
stituted by starbursts triggered by interactions of galaxies, 
as described in Menci et al. (2004) but not included in the 



SAM adopted in this paper. These starbursts would speed 
up the formation of stars in massive galaxies preferentially 
at high z (where the density of galaxies is larger). Such star- 
bursts would affect mainly the massive galaxies (due to their 
larger cross section for interactions) and would hence con- 
stitute the counterpart of the spheroids assumed to form at 
high-redshift in the PLE model. 

Of particular interest are the data by Lanzetta et al. (2002, 
solid diamonds in Figure 3), who found a monotonically in- 
creasing behaviour up to redshift 10. These data take into 
account also surface brightness dimming effects, which are 
likely to be serious at high redshift and which have never 
been considered before by any other group. In their most 
extreme case, the observations are as high as the values pre- 
dicted by the PLE curve uncorrected for dust. If confirmed 
by other deep surveys, the data by Lanzetta et al. (2002) 
could represent the most direct evidence in favor of a peak 
of star formation at high redshift. If true, such a peak would 
be problematic to explain for both PLE and hierarchical sce- 
narios. However, it also worth stressing that among the three 
sets of data calculated by Lanzetta et al. (2002) the most 
favored one by the authors is represented by the solid dia- 
monds with dotted error bars, of which the point at redshift 
z > 4 is in very good agreement with the PLE predictions 
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Figure 4. Evolution of tiie total comoving number density of galaxies brighter than 25.5 at rest-frame 1500A between redshift 2 and 
4.9, as predicted by by the PLE model {solid line) and by the hierarchical SAM of galaxy formation by Menci et al. (2002) {dotted line), 
and as observed by Steidel et al.(1999, open circles), Pozzctti et al. (1998, open squares), Lanzetta et al. (1999, open triangles), Chen et 
al. (1998, stars). This compilation of data has been taken from Somerville et al. (2001). For the sake of consistency with the data, in this 
case we assume a ACDM cosmology with Ho = 70km s^^ Mpc'^^. The theoretical predictions have been corrected for dust extinction. 
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Figure 5. Predicted fraction of the total present-day stellar mass as a function of redshift. Solid line: PLE model. Dotted line: SAM by 
Menci et al. (2002). 
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Figure 6. Upper panel: evolution of the stellar mass density as predicted by our PLE models (solid line), by the hierarchical SAM of 
galaxy formation by Menci et al. (2002) and as observed by various authors: Dickinson et al. (2003, square boxes), Cole et al. (2001, 
open circle), Brinchmann & Ellis (2000, open triangles), Cohen (2002, open squares), Fontana et al. (2003a, solid squares), Rudnick et 
al. (2003, open diamonds), Fontana et al. (2004, open pentagons). No mass cut has been applied to the predicted values. Lower panel: 
predicted evolution of the stellar mass density according to the PLE (solid curves) and SAM (dotted curves) by considering all the 
stars in galaxies with masses above three mass-cuts, namely M > IO^'^'^Mq (thick green lines), M > 10^'''^Mq (thick red lines) and 
M > IO'^'^'^Mq (thick black lines). The predictions are compared with observational values obtained by Glazebrook et al. (2004, solid 
circles) with the same criteria, i.e. by applying the three same mass cuts to the data sample. The values by Glazebrook et al. (2004) 
corresponding to the three cuts are plotted with the same colors as used for the theoretical predictions. 



but discordant with the SAM predictions. 
Also in the case of the high-redshift UV LD, the PLE and 
the hierarchical model used in the present work produce very 
different predictions, and the observations clearly allow us 
to discriminate between the two. 

Different indications seem to come from the UV luminosity 
density at z < 1. The prediction from the SAM by Menci 
et aL (2002) can nicely reproduce the data, whereas the 
PLE prediction is lower than the observations. At 2 = 0.2, 
where the lowest redshift observations have been performed 
by Treyer et al. (1998) at A = 2000A, the PLE models un- 
derestimates the data by a factor of ~ 2.2, whereas the data 
by Lanzetta et al. (2002) at 2: = 0.25 are underestimated by 
a factor of ~ 1.2. The explanation of this discrepancy is in 
part related to the fact that in the morphological classifica- 
tion of the PLE scenario we do not take into account nearby 
starburst galaxies, which can contribute up to the ~ 20% 
of the global star formation in the local universe (Brinch- 
mann et al. 2003). This would be enough to account for the 
discrepancy between the PLE predictions and the data by 
Lanzetta et al. (2002), but not for the data by Treyer et al. 
(1998). 

However, beside the missing contribution by starbursts, also 
the uncertainty in the B band LF normalization plays an im- 
portant role. The local B band LD adopted here for the PLE 



model is the one measured by Marzke et al. (1998), whose 
normalization is the lowest among the values provided by 
the most popular surveys (see Cross et al. 2001) and whose 
uncertainty could reach also factors of ~ 2. This fact could 
lead to a slight underestimation of all the LD values pre- 
dicted by the PLE model. 

The lower right panel of Figure 3 shows the observed evo- 
lution of the B band luminosity density compared with the 
predictions corrected for dust -extinction. At 2 < 2 the PLE 
and SAM curves are overlapping and both are in excellent 
agreement with the observations. At 2 > 2, the only avail- 
able measures are the ones by Dickinson et al. (2003) and 
by Rudnick et al. (2003) , none of which are accurately re- 
produced by any of these scenarios. In this case, however, 
the discrepancy is more critical for the PLE model than 
for the SAM. It is worth to stress that the combination of 
small field, cosmic variance effects, dust extinction and in- 
completeness are a non-negligible source of uncertainty in 
the data. Indeed, some of these effects cause also an under- 
abundance of massive galaxies as obtained by Dickinson et 
al. (2003) and a consequent underestimation of the stellar 
mass density with respect to the estimates by other authors 
(Fontana et al. 2003a, see also section 4.4). Also in the B 
band, absorption by dust significantly reduces the peak at 
2^4 — 5 due to ellipticals, although to a minor extent than 
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in the case of the UV band. In particular, the PLE model 
predicts a very narrow peak between redshift 5 and 4.8, 
corresponding to a time interval of ~ 60Myr. During this 
interval, the gas in spheroids is experiencing strong metal 
enrichment, consequently its optical depth is progressively 
rising to its maximum (see eq. 6) and the B band LD to 
its minimum. The fact that the peak is so narrow is due to 
the assumption that all spheroids start forming stars at the 
same redshift {zf — 5) and the star formation is completed 
after t ~ 0.3 Gyr. In a more realistic picture, the first galax- 
ies started forming stars before redshift 5 (see Giavalisco et 
al. 2004) and on a finite redshift range, so that the very nar- 
row peak would become larger and lower. Objects at high 
redshift which could be associated to a tail in the formation 
of galactic spheroids are the Lyman-break galaxies, which 
are usually detected at z > 3 and which show a large range 
of stellar population ages (Papovich et al. 2001, Shapley et 
al. 2001). In our picture, these galaxies can be associated to 
forming spheroids (see Matteucci & Pipino 2002), with to- 
tal stellar masses of the order of the Galactic Bulge. Other 
interesting objects are the submillimeter-bright galaxies, de- 
tected at z~ 2 — 3 and characterized by star formation rates 
of the order of 100- 1000 Ms/yr (Small et al. 2004). These 
galaxies have typical space densities of ~ 10~*Mpc~^, i.e. 
comparable to L* ellipticals (Blain et al. 2004). They ap- 
pear as massive as the largest spheroids observed locally 
and gas-rich (see Neri et al. 2003), and in the PLE picture 
they can be associated to a tail in the formation of massive 
spheroids. In a ACDM cosmology, the time lag between red- 
shift 2 and 5 corresponds to ~ 2.3 Gyr. This time-spread 
is consistent with what suggested by Bower et al. (1992), 
who found that in galaxy clusters the redshift range inter- 
ested by major spheroid formation could correspond to an 
age spread of A/orm ~ 2 Gyr. In the field, Bernardi et al. 
(1998) found a slightly larger age spread for large spheroids. 



i.e. A 



form 



3 Gyr. 



Another peak is predicted by the PLE curve a,t z ^ 
3.5, once the interstellar gas has completely been ejected by 
spheroids into the IGM, making the emission by the stars 
totally visible. 

Further observations in the B band at redshifts of 2-3 and 
beyond, within the reach of next generation deep galaxy 
surveys, could constitute a stringent test for PLE models. If 
the behaviour shown by present data should be confirmed 
by future surveys, this could constitute a strong evidence for 
galaxy density evolution, the process not taken into account 
in PLE models. 



4.3 The comoving galaxy number density 

In Figure 4 we plot the redshift evolution of the number 
density of bright galaxies. Such quantity is obtained by in- 
tegrating the rest-frame luminosity function at 1500 A, con- 
sidering only the objects brighter than the apparent mag- 
nitude limit of m,i5oo ~ 25.5. We consider only the redshift 
range between 2 = 2 and z — 5, i.e. the interval where the 
predictions provided by the PLE and hierarchical scenarios 
differ most. The observational data belong to various au- 
thors (see caption to Fig. 4 for further details) and have been 
all taken from Somerville et al. (2001). The observations in- 
dicate that most of the galaxy number evolution occurs in 



this redshift range: the number of bright galaxies is increas- 
ing by a factor of '^ 6 between redshift 2 = 5 and z ^ 2.8. 
The theoretical curves plotted in Figure 4 take into account 
dust corrections and represent the predictions according to 
the PLE (solid line) and hierarchical (dotted line) scenarios. 
The comparison between the theoretical predictions and the 
observations considered in this case indicates that the PLE 
scenario is inadequate to describe the number evolution of 
bright UV galaxies, since it systematically overestimates the 
observed number at all redshifts. We note that the disagree- 
ment between the PLE curve and the data is maximum at 
redshift z ~ 4, where the discrepancy is of a factor of ^^^ 5. 
On the other hand, the hierarchical scenario described by 
the SAM allows us to reproduce the observed trend with 
very good accuracy. It is worth noting that the study of the 
number density of bright UV galaxies represents an inter- 
esting test for the evolution of star forming galaxies at high 
redshift but, as well as the UV luminosity density, it does 
not provide any information about the formation of mas- 
sive spheroids, which most likely occurs in dust-enshrouded 
environments and are thus invisible in the rest-frame UV. 
Furthermore, if at redshift 3-4 there was already a significant 
number of massive galaxies containing old stars, generating 
red spectra, such population would be certainly missed by 
UV galaxy surveys. A fruitful test for the identification of 
the number of massive galaxies at high redshift is the study 
of the evolution of the stellar mass density. 

4.4 The evolution of the stellar mass density 

Figure 5 shows the redshift evolution of the stellar mass 
fraction as predicted by the PLE model (solid line) and by 
the SAM (dashed line). Each curve is normalized to the 
value for the stellar mass density predicted at the present- 
day. This figure is helpful to understand what percentage 
of the present-day stellar mass is in place at any given red- 
shift according to the predictions of the two scenarios. The 
two curves have a very different behaviour: according to the 
PLE model, nearly half of the stars observable today are 
already in place at z ~ 4, corresponding to 1.63 Gyr after 
the big bang for the cosmology adopted here. This is due to 
the stellar mass produced in spheroids. The increase from 
2 = 4 to 2 = is due to quiescent star formation in spirals 
(CM03). At 2 = 1, corresponding to an age of the universe 
of 6.2 Gyr, the PLE model predicts that 85% of the present 
stellar mass is already in place. 

According to the hierarchical SAM, the buildup of the stel- 
lar mass occurs progressively, with half of the total stellar 
mass in place at 2 ~ 1.2, i.e. 5.42 Gyr after the big bang. 
By 2 ~ 1, the SAM predicts that nearly 60% of the total 
present stellar mass is present. 

Figure 6 shows a comparison between the stellar mass den- 
sity as observed by various authors and as predicted by PLE 
models and by the SAM. This comparison demonstrates 
that, owing to the extreme differences between the PLE and 
SAM predictions, the observation of the stellar mass density 
constitutes another very helpful strategy to distinguish be- 
tween the PLE and the hierarchical scenario. 
In the upper panel of figure 6, we show the evolution of the 
stellar mass density by considering galaxies of all masses, 
namely no mass cut has been applied to the predicted val- 
ues. The theoretical predictions are compared with obser- 
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vational estimates by various authors (for further details, 
see caption of Fig. 6) In general, the main sources of uncer- 
tainties in the data are dust extinction and cosmic variance 
effects due to small field. The data by Fontana et al. (2003a) 
are taken from a large volume and are corrected for dust ex- 
tinction. However, as emphasized by the authors, they may 
still suffer for incompleteness on the bright end of the mass 
function. To estimate to what extent these effects could alter 
the real values is difficult: for instance, the amounts of dust 
can vary considerably from one galaxy to another. Also the 
cosmic variance effects are in principle difficult to evaluate. 
It is worth noting that all these effects conspire to lower the 
estimates of the stellar mass at redshifts larger than 1: for 
these reasons, it is safe to consider the data as lower limits. 
The PLE and SAM curves are both in reasonable agree- 
ment with the data within redshift z < 1.5. At redshifts 
higher than 1.5, if we consider the predicted total stellar 
mass the PLE model presents a noticeable discrepancy with 
the observations: if we consider the central values estimated 
by Fontana et al. (2003), the discrepancies between obser- 
vations and PLE predictions are by factors of 3 — 6. On the 
other hand, on average, the SAM predictions seem to show 
a good agreement with the observed values. 
In the lower panel of figure 6, we show the predicted evolu- 
tion of the stellar mass density according to the PLE (solid 
curves) and SAM (dotted curves) and by considering all the 
stars in galaxies with masses above three mass-cuts, namely 
M > IO^^-^Mq (thick green fines), M > IO^^'^Mq (thick 
red lines) and M > lO^^'^M© (thick black lines). Such pre- 
dictions are compared with observational values obtained by 
Glazebrook et al. (2004) with the same criteria, i.e. by apply- 
ing the same three mass cuts to the data sample. The values 
by Glazebrook et al. (2004) , corresponding to the three cuts, 
are plotted with the same color as used for the theoretical 
predictions. The adoption of the mass cuts is very helpful 
in establishing a full correspondence between observations 
and theoretical predictions, and to have a very clear picture 
of the number of massive galaxies that the PLE and hier- 
archical scenarios predict at any redshift, respectively. If we 
compare the PLE predictions with the data calculated with 
the three cuts, we notice that the agreement between data 
and predictions does not improve and that the PLE model 
in general tends to overestimate the stellar mass density in 
massive galaxies, in particular at redshifts z > 1. 
If we compare the SAM predictions to the data, we no- 
tice that the hierarchical picture can reproduce the observed 
data with the three cuts up to redshift z ~ 1.2, whereas at 
higher redshift it tends to underestimate the observations. 
The disagreement is particularly strong for the highest mass 
cut (M > 1O"'-*M0). This shows that at redshifts z > I, 
according to the SAM predictions, the bulk of the stellar 
mass resides in objects with masses M < 10^"'^ M0. These 
small objects would be too faint to be visible by any current 
high-redshift survey. Also in this case, this problem is allevi- 
ated by considering the effect of interaction-driven starburst 
in massive galaxies at high-redshift, (see Menci et al. 2004), 
which would increase the fraction of stellar mass already in 
place at z = 2 to a value around 0.3 of the present mass 
density. 

It is very interesting to see how, by means of ACDM cosmo- 
logical numerical simulations, Nagamine et al. (2004) find 
a strong discrepancy between the predicted and observed 



amount of stellar mass at redshift z > 1.5. Their simula- 
tions indicate an excess of stellar mass with respect to ob- 
servational estimates at high redshift, in analogy with the 
result of the PLE model considered in this work. This is an- 
other indication suggesting that the global star formation of 
the universe may have proceeded in the past at levels some- 
what higher than predicted by semi-analytical models, and 
it confirms that effects such as dust obscuration and cosmic 
variance may still seriously prevent us from having a clear 
picture of galaxy evolution at redshifts z > 1. 
Recently, the Great Observatories Origins Deep Survey has 
provided evidence for a population of galaxies showing dis- 
torted morphologies and with ongoing merger activity lo- 
cated at z > 1.5 (Somerville et al. 2004). The number den- 
sity of such bright objects is underestimated by current hier- 
archical SAMs and overestimated by PLE models. To assess 
the role of such galaxies in the stellar and metal budget 
would be of primary interest in order to have further crucial 
hints on the evolution of galaxies at redshifts larger than 1. 



5 CONCLUSIONS 

In this paper we have studied the evolution of the cosmic 
star formation, the galaxy luminosity density and the stel- 
lar mass density by means of two opposite galaxy evolution 
pictures: the pure-luminosity evolution model developed by 
CMOS and the semi-analytical model of hierarchical galaxy 
formation by Menci et al. (2002). The former predicts a 
peak at redshift z = 5, due to intense star formation in 
ellipticals, followed by a phase of quiescent and continuous 
star formation occurring in spiral galaxies. The SAM pre- 
dicts a smoother behavior, following the gradual build up 
of galactic DM halos through repeated merging events. The 
aim was to derive constraints on the relative importance of 
different physical processes - like the dependence on mor- 
phology of the star formation history, the density evolution 
of the galaxy population, the impulsive star bursts - in de- 
termining the observed properties of the galaxies. 

We have shown that the evolution of the cosmic star 
formation rate density in the two models behaves quite dif- 
ferently. However, the integral of the cosmic star formation 
rate at redshift z < 1, probed by the stellar mass density 
evolution in this redshift range, are in good agreement. This 
ensures that the total amount of stars formed along the star 
formation histories are similar (and in agreement with the 
observations) . To probe the rate of star formation at differ- 
ent cosmic epochs we investigated the luminosity density in 
the UV and B bands, where the emission is dominated by 
young, short-living massive stars. The comparison with the 
available data shows that: 

1) At redshift z > 4, the SAM tends to underesti- 
mate the observed UV luminosity density which, as several 
current surveys indicate, is a non-decreasing function of z. 
On the other hand, the PLE predictions can fairly account 
for such observed trend. If future surveys will confirm such 
behaviour, this could indicate that some fundamental pro- 
cesses should be inserted into SAM to boost the star forma- 
tion at high redshifts. An example of such a process could 
be the interaction-driven starbursts suggested by Menci et 
al. (2004). 

2) In the B band the PLE model tends to overestimate 
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the observed luminosity density at 2 > 2.5 by a factor in- 
creasing with z. This is the consequence of placing a rapid 
formation of all the elliptical galaxies at z « 5. While dust 
extinction and incompleteness severely affect the compari- 
son with present data, if future observations will not indi- 
cate a substantial growth of the B-band luminosity density 
for z > 2.5 — 3, this would point toward a galaxy density 
evolution, the main process not included in PLE models. 

3) At low redshift {z < 1), the local UV luminosity 
density predicted by the SAM is about 2 times larger than 
those arising from PLE models. This is because in hierar- 
chical scenarios at low redshift the small mass galaxies still 
retain a significant fraction of their cold gas reservoirs, while 
the massive ones have already exhausted most of their fuel at 
high redshift, since the latter are formed from clumps orig- 
inated in biased high density regions of the cosmic density 
field. In hierarchical models, at low z the contribution of low- 
mass galaxies sustains the global star formation rate above 
the value obtained in the continuous, passive evolution PLE 
models. The above discussion shows that, while the local 
J and K observations will hardly contribute to discriminate 
between the two scenarios, accurate measurements of the lo- 
cal UV luminosity density would be effective in constraining 
the models. 

4) The observed evolution of the comoving number den- 
sity of bright galaxies at redshift 2 > 3 is well reproduced 
by the hierarchical SAM, whereas, for the set of data con- 
sidered here, the PLE overestimates the observed densities 
by factors between 2 and 5. 

5) The stellar mass density constitutes a complemen- 
tary probe for the PLE and hierarchical scenarios. In gen- 
eral, both the PLE and hierarchical predictions allow us to 
reproduce the observed stellar mass density evolution up to 
2 = 1. At 2 > 1, the predicted stellar mass densities di- 
verge, with the PLE predictions remaining almost constant 
up to redshift 2 ~ 4 and the SAM predictions continuously 
dropping with increasing 2. Without any mass-cut on the 
theoretical predictions, the PLE model overestimates the 
data by factors of 3-6. If we calculate the stellar mass den- 
sity evolution and apply the three mass cuts, as performed 
by Glazebrook et al. (2004), in general the discrepancies be- 
tween the PLE model and the observations at 2: > 1 do 
not reduce. On the other hand, the hierarchical picture un- 
derestimates the observations for all the three values of the 
mass-cuts at redshifts 2 > 1.2. This is related to the fact 
that, at redshifts 2 > 1, according to the SAM predictions 
the bulk of the stellar mass resides in objects with masses 
M < 10^°'^ M0. These small objects would be too faint to 
be visible by any current high-redshift survey. Also in this 
case, the discrepancy between the hierarchical model and 
observations is partially alleviated by introducing a popula- 
tion of high-redshift starbursts in massive galaxies (Menci 
et al. 2004) , which would bring the mass density at 2 = 2 to 
values around 1/3 of the local value, in much better agree- 
ment with the data but still well below the PLE predictions. 
Thus, in principle, more precise observations of the stellar 
mass density at 2 > 2 will be able to discriminate between 
the PLE models and the SAM including starbursts at high 
2. On the other hand, some indications against hierarchi- 
cal formation of elliptical galaxies is provided by chemical 
constraints, in particular the increase of the [Mg/Fe] ratio 
with galaxy luminosity (Pipino & Matteucci 2004, Thomas 



1999 ). This fact indicates that the most massive ellipticals 
stopped forming stars before the less massive ones. All of 
these facts together will have to be taken into account even- 
tually before drawing firm conclusions. 

As forthcoming work, to investigate star and massive galaxy 
formation at high redshift we will use other diagnostics, such 
as IR and submm emission. 
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